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I. INTRODUCTION AND BACKGROUND 

Previous studies ^ ^ have identified 

(1) A range of central scientific questions concerning combustible 
systems— which are best studied with the aid of a Space Shuttle 
laboratory facility. 

(2) "Flame propagation and extinction in clouds of porous particulates" 
as one of a select group of important combustion experiments for 
such Space Lab investigations. 

(3) An experimental approach, apparatus, and procedures for the study 
of item (2) , above 

(4) A theoretical approach to the necessary understanding of the 
proposed experiments of item (2), above. 

The current research effort is designed to advance the theoretical structures 
needed for the predictive analyses and interpretations for "flame propagation 
and extinction for clouds of porous particulates." An effort has also been 
made to further advance related combustion theory of significance to reduced 
gravitational studies of combustible media. 

Theoretical studies - ^ of flame propagation through clouds of porous 
particulates are generally couched in terms of sets of conservation equations 
and appropriate boundary conditions. Central to the requirement that a flame 
propagation theory also be useful for flame extinction analyses, is the re- 
quirement that the boundary conditions be nonadiabatic. Nonadiabatic boundaries 
are required for both autoignition theory and for extinction theory. Processes 
that have been considered include, among others, pyrolysis and vaporization of 
particulates, heterogeneous and homogeneous chemical kinetics, molecular trans- 
port of heat and mass, radiative coupling of the medium to its environment. 
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and radiative coupling anong particles and volume elements of the combustible 
medium. 

In general, the needed theoretical structures provided thereby should be 
capable of quantitatively representing: 

(a) Steady state flame propagation 

(b) Flame extinction limits 

(c) Small particle flame behavior and limits 

(d) Large particle flame behavior and limits 

(e) Effects of apparatus scale 

(f) Effects of system transport properties 

(g) Effects of variations in kinetic characteristics of combustibles 

(h) Autoignition conditions . 

It is expected that the results of this investigation will help to provide 
an analytical model that can be used to 

(a) Aid in experiment selection and design 

(b) Understand in-space: experimental results 

(c) Provide a predictive approach to further understanding of combustion 
for g >_ 0. 

In addition to the above cited studies, we have performed related theo- 
retical examinations of selected areas in the field of combustion where reduced 
gravity studies appear promising and where investigations to date appear inade- 
quate. These latter analyses attempt to identify the critical aspects of the 
associated combustion processes that appear to require investigation under 
zero-(g) conditions. It is anticipated that these ab initio theoretical studies 

will be useful to our understanding of fundamental phenomena associated with 

(7 8 ) 

oscillatory combustion. 
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II. GENERAL THEORETICAL CONSIDERATIONS 


For a cloud of particles at reduced gravitational conditions, we consider 
a flame to be propagating in the negative x-direction. The one-dimensionalized, 
nonadiabatic energy conservation equation can be written for the gas phase 
components of the system 


( 1 ) 


-f- [G(c T + h ] 
dx v g g v J 


d 2 T 


\ 


dx 



+ 


q + CS p q p N) - Z L ± 


where 


x - axial direction 

Cg - specific heat of the gaseous medium 
G *!■ mass flow rate 

h v - heat of vaporization of the condensed phase 
Tg - Gas temperature 
A - molecular conductivity 

(dq /dx) - axial component of the divergence of the radiative flux density 
K 

q' 1 ' - heat release rate associated with gas phase chemical reactions 

q^' - local rate of a particle's rate of heat rele&se (or heat 
absorption) 

Op - local concentration of particles 
EL^ - Sum of the radial heat loss rates due to all transport processes 


Equation (1) provides for either/both exothermic/endothermic processes 
within a particles spatial regime, as well as for gas phase chemical kinetics. 

The two-dimensionality and further complexities of the radiative-conduc- 
tive coupling which may be important have been taken to be separable. 

This latter simplification may not always be useful and an extensive examina- 
tion of this issue is provided in a later section of this report. 
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Gasification due to individual particles may be accounted for by the 
expression 


(2) m - A expC-E /RT ) 

p s r v p 

Within a' propagating flame's structure, the total gaseous mass flow rate, G, 
is not constant, but increases due to the net effect of all gasification 
processes. Accordingly, 

«> Jf - V*” ‘-V RT p)n p 

provides the matching condition which couples the vaporization reaction rate 
of the particle cloud to the gaseous mass flux increase. The total mass flux 
rate, G c> is given by 



In previous discussions, we have noted the general approaches to numerical 
solution of a system of equations such as Cl) — (5) subject to appropriate boun- 
dary conditions. 

Several interesting conditions limit the ranges of parametric conditions 
within which quasi-steady flame propagation can be achieved. The first of 
these pertains to an upper limit on the ambient temperature conditions. 

At sufficiently high ambient temperature conditions, autoignition occurs. At 
or above the autoignition condition, no quasi-steady flame propagation is pos- 
sible (the system is supercritical) . Other parametric conditions which may 


limit existence of quasi-steady flame propagation derive from either energy 
loss to the environment or from inadequate coupling of the two-phase heat and 
mass transfer mechanisms which are implicit in equations (1) — (5) . Both auto 
ignition and flame extinction limits are discussed in the following portions 
of this, report. • 
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III. AUTOIGNITION 

Ignition and "explosion" supported by a uniform gaseous system has been 
treated theoretically and reviewed previously . The elements of such a. 

'theoretical structure already exist in appropriate nonadiabatic gas phase flame 
propagatipn theory. One has only to set all convective flows to zero and 

examine the classes of solutions for "nonflowing gaseous systems". This has 
been done. In a similar way, one may reexamine the limiting case for auto- 
ignition for a uniform dust cloud of combustible particulates. Data accumulated 
by the U. S. Bureau of Mines v ' provide a range of empirical "explosibility 
indices" for clouds of particulates. Autoignition temperatures which are theo- 
retically comparable to those for premixed gaseous systems are derivable from 

(13) 

Cl) — (5), and related considerations. 

The theory of the autoignition of single particles of metal and of coal 

has received considerable attention. ^ However, experiments show that 

the cooperative effect of the particles in a cloud has a significant influence 

on the ignition temperature and that the ignition temperature can also be 

dependent on the system’s physical^^’^^ parameters. Meese and Skifstad have 

(19) 

treated some features of the ignition of clouds of boron particles. Their 

definition of ignition temperature is rather specific to the type of metal 
under consideration (boron) and computations are carried out for ignition times. 
The definitions and the aims, as well as the mathematical treatment developed 
here, are different and more general. Rumanov and Khaikin developed an ap- 
proach^^ to ignition of particle clouds. The present work uses methods of 
the phase-plane^'*'^ to describe ignition conditions resulting from equations 
(1) — (5), with no convective flows. 
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We consider the problem of autoignition of a cloud of solid particles 
uniformly dispersed in a gas with one component of which it can react exother- 
mally. The heat generated by the gas-solid reaction heats the particle and 
the surrounding gas is heated by particle-gas heat transfer. The mass of gas 

. in the cloud loses heat to: the ambient at the cloud boundary (or a vessel 

( 12 ) 

boundary as in laboratory experiments ). The particles can also directly 
exchange heat with the ambient through radiation; radiative transfer is inclu- 
ded only to indicate its importance in some cases, but is not treated in 
detail here. 

We assume the following: 

1) The system is "homogeneous" in the sense of having uniform properties 
across the cloud. (This is similar to the classical Semenov theory of auto- 
ignition in gaseous systems.) We adopt a two-temperature model with one tem- 
perature for the particles aid another for the gas in the cloud. 

2) A heterogeneous reaction occurs at the surface of the particles (for 
example, if the cloud consists of metal particles in air, this could be a sur- 
face oxidation process) . We can generalize the assumed surface reaction to the 
case of vaporizing particles in subsequent analyses. For the case of vaporizing 
particles, gas phase autoignition conditions similar to those previously 
discussed are obtainable. 

3) The induction period is sufficiently small so that negligible con- 
sumption of fuel or oxidant occurs. 

4) Heat transfer is represented by Nusselt's number type coefficients. 
Radiation transfer is negligible. 

5) The effects of natural convection are negligible. (These effects can 
be important for large temperature differences and cloud sizes.) 

6) All physical properties are constant. 
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The "homogeneous" model permits us to write two energy equations for the 
cloud, one for the particle, and the other for the mass of gas in the cloud. 
They are 


( 6 ) 

and 

( 7 ) 


dT 

M C -rr 2 - 
p p dt 


dT 

M C -r-£ 
g g dt 


Q - a. S (T - T ) - L 
T 1 P P g i r 


Na,S (T - T ) - a.S(T 
1 P P g 2 g 


U 


where 


M - Particle mass 
P 

0 ^ - specific heat of the particulate material 
T^ - ambient temperature at the boundaries 
T p - temperature of a particle 
t - time 

Qp - net rate of heat generation per particle 

- coefficient for heat exchange between particle and gas 
a ^ - coefficient for; heat exchange between cloud and surroundings 
L r - total rate of heat loss by radiation per particle 


N - total number of; particles 

S - internal surface area of the vessel containing the particulate 
cloud 

S - particle-surface area 
P 

Yq - mass fraction of reacting gas component. 


The nature of the above equations suggests a phase-plane representation^'*'^ 
The advantage of such a method is that the local stability of a steady state 
can be inferred from the classification of the corresponding singular point in 
the phase-plane, even where the defining equations are nonlinear. 
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In the T - T phase-plane, the singular points (and the corresponding 
PS 

steady states for the system) are given by solving the following equations 
simultaneously: 


(8) NctiSpCTp - T g ) - a 2 S(T g - Tj = 0 


Q - a.S (T - T ) - L = 0. 
P 1 P P g r 


Simplifying, 


% = (T P " T »> + L r 


a, S 

— i_2 (-p 

Na S p 

L. IP 


TJ + V 


If Qp is an Arrhenius-type rate with an exponential dependence on t , 
the particle temperature, then Eq. (10) would give the familiar Semenov-type 
distribution of singular points (steady states). The nature (or local stability) 
of each of these steady states is determined by the first partial derivatives 
with respect to the phase-plane variables of the expressions on the right-hand 
side of equations (6) and (7) evaluated at the singular point . For an 
Arrhenius -like Q^, it can be shown, by the methods set out in (11) that the 
two lower singular points would be a node and .a saddle point. The critical 
condition for ignition occurs when the two points merge and is given by 


a.S 9L 

l p + — £ 

Na. S 3T 

+ ~kL P 

a l s 
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The simultaneous solution of Eqs. (10) and (11) (if is known as a 
function of Tp) will give the autoignition condition, that is, the value of 
T^ for autoignition. Along with Eq. (9), the following expressions, are obtained 


.for. T and T : 
00 g 


T m = T - 

00 P 


«p ~ V M( % - L r> 

Vp ' “2 S 


(13> T g ’ T p " ( VVp > + V a l V- 


The examination of these equations requires specification of the forms 


for 0^ and L r . The effects of radiative loss would be important in high- 
temperature situations and the expression for L r should be dependent on the 
cloud density. However, for the present L r is set equal to zero. 

Some general results can be obtained by examining limiting cases with 
the help of Eqs. (11) — (13.) . We assume that Qp can be described as a simple- 
surface reaction by: 

(14) Q p = S p Ah c pY 0 exp (-E/RTp) . 

3/2 3/2 

Cloud-particle density is proportional to NSp /S if we assume spherical 
particle and cloud geometry. 

We may now consider several alternate cases : 

Case (a) : When No^Sp/c^S 1, Eqs. (11) and (12) reduce to 

U5 > iVexpc-E/i " sr -Hprs- 

p p 0 c 


s 3/2 v 


2 |»s 3/2 Jl S 1/2 j [“P Y 0 h cJ 
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: and 

(16) T^ = T (1 - RTp/E) . 

From Eqs. (15) and (16), we can deduce that for the same particle and 
cloud size, the temperatures at ignition (T^ and T^) decrease with increase in 
cloud density. For the same density, the ignition temperatures decrease with a 
decrease in particle size and with an increase of cloud size. 

Case (b) : When Not^Sp/o^SCk 1, equations (11) and (12) reduce to: 

(17) e^C-E/RX ) - 

p 0 c 

(18) T - T (1 - RT /E). 

P P 

Equations (17) and (18) represent the single-particle limit in that the 
cloud parameters are absent from the ignition conditions. The particle-size 
dependence comes in only through such dependence of a^, the particle-gas heat 
transfer coefficient. For. stagnant conditions, is inversely proportional 
to the particle diameter and the ignition temperatures would decrease as parti- 
cle size increased. This also suggests that ignition is likely to be observed 
as occurrences on individual particles rather than as cloud phenomenon. 

The limiting-case results above are qualitatively supported by some of 
the experimental results of Cassel and Liebman^"^ on ignition of metal dust 
clouds. Their Figs. 2 and 3 show that the ignition temperatures fall with 
increase in dust concentration for a given particle size. The different 
dependence on particle size in the two limiting cases is also seen by the cross- 
over of their ignition temperature-cloud density curves for different particle 
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In Case (a), the ignition temperature at constant cloud density is 
independent of particle size, While a decrease with increase of particle size 
is shown for Case (b) . 

From this brief examination, it appears that the phase-plane technique can 
be fruitfully employed in determining autoignition conditions for dust clouds. 

The ignition conditions derived show qualitative agreement with certain experi- 
mental measurements on metal-particle clouds. The method can be extended to 
study systems with different chemistry and including radiative cooling. However, 
it appears necessary tocarry out autoignition experiments at reduced gravatational 
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conditions if energy transport is to be restricted to conduction and radiation 
and if. truly homogeneous clouds of particulates are to be studied. 
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IV. FLAME EXTINCTION 

In the study of premixed, gaseous flames it is found that extinction 
limits for flame propagation are associated with a number of critical condi- 
tions beyond which integrated heat release rates (for the flame structure) 

cannot match the integrated heat loss rates (integrated over all boundaries) 

(21-23) 

sustained by the flame structure. 

Equations (1)— (5) characterize a gas phase medium which, if it is to 
support flame propagation through. a cloud of porous, vaporizable solids, must 
be supplied with adequate gaseous fuel for combustion. Accordingly, it is 
immediately clear that the. boundaries of interaction for the gaseous medium 
(of this two-phase system) are important in ways which don't apply for pure 
gas phase media. These include: 

(1) Heat losses from the gaseous medium to the particulates are required 
to support the vaporization process which is necessary for provision 
of gaseous fuel. 

(2) The rate of heat loss sustained by the gaseous component (of the 
two-phase system) includes, as boundary losses, both the apparatus 
boundaries which enclose the particulate cloud and the sum of energy 
losses to the vaporizing particulates. 

(3) Hot, optically inactive gaseous components of flame gases may 
indirectly sustain local flame radiation loss rates via the radia- 
tive coupling among particulates as well as the radiative interaction 
between the particulate cloud and the cold boundaries of the reactive 
medium. 

(4) yaporizahle materials for which (E^) is very high and for which par- 
ticle sizes are large (at constant total particle mass concentra- 
tions) will generally fail to support flame propagation. For such 

a system 
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(a) volumetric kinetic vaporization rates are very small, except 
at the highest temperatures (Eq. (2)). 

(b) volumetric heat transfer (and vaporization) rates from gas to 
particles at the highest temperatures are still very slow, due 

; ' to the small number of particles and the fact that these rates 

depend upon (T - T ) values; 

S P 

(c) For a broad range of "prevaporization temperatures", the parti- 
culate cloud sustains radiative losses to the boundaries, yet 
contributes little gaseous fuel to the gas-phase flame system. 

Ihese observations are generally in accordance with the qualitative 
features observed' for flame propagation limits. However, quantitative appli- 
cation of the approach discussed herein requires observation of flame propaga- 

(24) 

tion and extinction data obtained under conditions of weightlessness . We 

(24) 

have noted previously that flame propagation data obtained at normal 
gravitational conditions suffer from a number of necessary deficiencies. 
Additionally , when flame propagation rates are at all comparable with particle 
settling rates, apparent flame extinction (at normal gravitational conditions) 
cannot be attributed primarily to the transport and kinetic processes discussed 
here. Accordingly , the limiting two-phase coal flame propagation conditions 
reported in reference (6) need not at all correspond to what we may observe 
under zero gravitational conditions. Preliminary observations with lycopodium 
powder-air flames (at reduced gravitational conditions) are in conformity with 


this observation. 



The Conductive-Radiative Mechanisms 


of Dust Flames 

1. Introduction 

Previous theoretical models for propagating coal-dust flames have 

considered transport mechanisms based upon radiative heat transfer from 

16 25 26 

the hot flame region to the incoming coal /air mixture. * * These theories 

appear to be applicable only to devolatilized char flames wherein the heat 
release rates are slow. In general, the rapid volatilization in coal-dust 
flames leads to the formation of thin flame reaction zones where temperature 
rises steeply. The existence of such thin flame fronts (y 1 cm) suggests 
"the importance of gas phase diffusional (conduction) processes in such 
flames".*’ In the early theoretical model of Smoot, Horton, and Williams** 
radiative processes are ignored. More recent studies do not consider the 
general conductive-radiative problem for the reaction zone. 

If the optical thickness of a flame front based on the inverse of the 
Rosseland mean free path is either very large or very small, the effect of 
radiative processes on the structure of flame reaction zone is expected to 
be small. In the case of; a very large optical thickness, this result is 
due to the fact that the radiative flux rate is small. In the case of very 
small optical thickness, this result is due to the fact that the divergence 
of the radiative flux rate is small. In transparent gaseous flames the 
Rosseland mean free path is infinite and the corresponding optical thick- 
ness is zero. The effect of radiation on such flame structure is therefore not 
important. In dust flames, the Rosseland mean free path is neither infinite 
nor zero. Thus, the roles played by radiative transport in dust flames 
requires further examination. 
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Methods of aialysis of radiative transfer problems vary greatly from 
medium to medium, depending upon the different optical properties involved. 
It is therefore necessary to briefly review the optical properties of a 
typical dust medium. This will be useful in the study of the conductive- 
radiative theory of dust flames to be formulated. It will also serve to 
guide us in selecting the ; analytical tools needed to formulate the theory. 

In previous experimental studies^, coal particles of sizes between 
3 y to 38 y were used. Due to its ready availability, we shall use the 
optical properties of carbon black particles (refractive index, 2.00—1.20 i) 
as an approximation to those of coal particles. 

Calculation of Mie solutions for spherical particles with refractive 

27 

index, 2.00—1.20 i, was carried out by Fahimian. Let the dimensionless 
particle size, x, be defined as, 

x = 2irr/X, . 


where r is the radius of the spherical particle and X is the wavelength 

of radiation. Let K , K , E be the extinction coefficient, the scattering 

C S d 

coefficient, and the absorption coefficient, respectively. The dimensionless 
form of the (extinction) coefficient is called the (extinction) efficiency, 

X, which is defined as, 



W 



= K /C A 
s n 


X 

a 


= K /C A 
a n 



where C n is the number density of particles and A is the geometric cross- 
section area of individual particle. The Mie solution predicts the relation 
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27 

uetween the efficiency and the dimensionless size. The results are repro- 
duced in Fig. 1 in the forms of - x and X /X - x curves. 

° t s a 

For dust media of kncvm mass concentration of solid particles, C , 

m 

and the size of particles, the spectral extinction coefficient, K^CA), and 
the spectral single-scattering albedo, tTq(A), can be readily calculated, 


K„(A) = X (x)C A = X (x) 


V X) “ K. " 1 X 

t s 


fC 1 4 


fl 


+ 1 


where p is the density of the solid phase and V is the volume of the solid 
particles. For spherical particles, 


A _ _3_ 
V 4r 


In the ranges of particle sizes and of wavelengths considered, the value of 
x. is greater than unity and the variation of -X is moderate. With decreasing 
particle size, the decrease in X is moderate and the increase in the 
area/ vo luma, ratio, is more pronounced . Consequently,, the extinction- coeff in - 
dent increases with decreasing particle size under constant mass concentra- 
tion. This property, vhieh is displayed in Table i and Fig. 2, may effect 
the dependence of the burning velocity of dust flames on the mass concentra- 
tion and the particle size (see Fig. 4). We shall come back to this point 
later in our discussion. 

The other significant (optical) property of dust media with particles 
is this size range and radiation field (in this wavelength range) derives . 
from the fact that the variation of the extinction coefficient and the 
.single-scattering albedo with wavelength is small. It is therefore possible 
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to consider the idealization of dust-containing media as gray media. A 
good estimate of the error involved in such a gray approximation is given by 
the difference of the Planck mean from the Rosseland mean extinction coeffi- 
cient, i.e., 



K 


tP K tR 


K 


tR 



Such an estimate is shown in Fig. 3 and supports fully the use of the gray 
approximation. ; 

Thus, the adoption of a mean extinction coefficient is made possible. 


3 . The conductive-radiative transport mechanism 

We now consider the relation between heat transfer processes and 
reaction zone structure and burning velocity. Differences in particle and 
gas temperature will be ignored. 

The energy balance may be written as 


(24) 


a 


d 2 T 

dx 2 


p a C 


V* (htF) + H - w 


dr 

dx 


tTt 7 u ’ 


where a is the thermal diffusivity, 

_ 00 

4ttF(= 4tt ^ F^dv) is the radiation flux, 

H is the conduction and radiation heat loss, 
p is the density of medium, 
w is the burning velocity, 

QC is the heat of reaction, 

1/t is the rate of reaction. 


For an adiabatic flame, H = 0. One may consider equation (24) the equation 
for the unknown T since F is a known functional of T. The equation is an 
integral-differential equation with integration over both frequency and 
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physical space. Such an equation is difficult to solve. However, once the 
flame is considered as a gray medium, it becomes possible to remove the 
need of integration not only over frequency space but also over physical 
space. The approximation which removes the need for integration over 
physical space is called the diffusion approximation. 

The monochromatic radiation flux, 4irF^, may be expressed approximately 
in terms of the monochromatic Planck's function, 'B , and the monochromatic 
mean intensity, J ., 


(25) 


4-tt, 


dJ . 


v 


9(1 "iV dTv 


l-TT, 


9(1-3 V 


dT v 



where is the optical depth and 7T^ is the asymetric factor of the scat- 
tering phase function. Equation (25), which was derived by using the 

op 

exponential kernel substitution method , reduces in the limit of radiative 
equilibrium in gray media to 



J = B 


V V 



In the limit of conservative scattering, 


*0 = 1 • 


and the classical Eddington approximation is obtained 


dJ 


(26) 


F = - 
V 


V 


T/i 1 . dx 

3 ( x ~ 3* F 1 ) V 


It has been shown that in the general case this approximation, equation (25) 
remains quite accurate while the Eddington approximation, equation (26) , fails 
in the case of 


TTq « 1 


J i b 
v v 


and 
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With K t (V) = K t ^(~K) , equation (25) can be integrated with respect to 
frequency and becomes 


(27) 


4-ir, 


F = - 


9(1 -y 7r x ) 


’ dJ 
dT 


1-TT, 


9d- y 


dB 

dt 


where B = /Bd = — T atl d a is the Stefan-Boltzmann constant. Radiation 
0 v ir 

flux is then expressed, according to the diffusion approximation in terms 
of temperature and mean intensity. One needs another differential equation 
which together with equation (24) describes the temperature distribution and 
the mean intensity distribution. This equation is the emission-absorption 
equation. 


(28) yy = (1-7T 0 ) (B-J) . 

Eliminating F among equations (24) , (27) , and (28) , we have the governing 
equation for T and J, 


(29) 


“ TT - iTc a-VKlB(T)-j) - » § + - 0 

dx 'a N 


(30) 




9(1 -y ^i) dt" 9(1 -y^) dT 


With additional equations to describe x(T), and with appropriate boundary 
conditions, these equations can be integrated to determine the eigenvalue, w. 


4. Burning velocity of coal-dust flames 

Without obtaining explicit solutions to equations (29) and (30) we may 
examine them to make a qualitative coraparision of a conductive-radiative 
theoretical model with the pure conductive theoretical model of coal-dust 
flames . 
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Experimental flame velocity measurements were reported by Horton and 
6 

Smoot. "In general, a lean flammability limit is evident around 100-200 mg 
coal dust per liter of air, while maximum flame velocities of 20—35 cm/sec 
at 300 to 600 mg/2 and rich limits well over 1000 mg/2 were also observed. 

In addition, an increased particle size markedly decreased the lean-side 
and peak flame velocities while causing the peak to occur in a richer flame. 
However, on the very rich side, a large size increased the flame velocity 
... . " A typical example of experimental results together with their theo-r 
retical results is reproduced in Fig. 4 . 

Consider the small particles of 10 )J at concentrations greater than 
500 mg/2. The optical thickness of the flame front is greater than unity. 
At 1200 mg/2 the front optical thickness is about 3, i.e., a moderately 
large : optical thickness. In the limit of large optical thickness, J - B 
becomes small, i.e., 

J = B 


and equations (29) and (30) reduce to 


(31) 


?a C 


- d 
dx 


dT , ' ' 16CTT dT 

dK 3(l-|TT 1 )K dx 


.w'f+m.-- o 

dx x(T; 


The effective "conductivity" is 


29 


k + 


16a T“ 


3(1 tt 1 )K 


The heat flux in a pure conductive model is 


k dx ’ 



while' in a conductive-radiative model the heat flux is 


“ k .( 1+ N- 1 dx » 


where — • = 
N 


16aT" 


3(l-y Tr 1 )kK 


is a dimensionless parameter measuring, the relative 


importance of radiation and conduction. Thus, consider the example where 

: T = 1200°K 

k = 0.0782 w/nK for air 

K(l-yTr^) = 300 m ^ (corresponding to roughly 10 y 
particles at 1200 mg /£.), 


we find 


- = 22 . 

N 


Thus the pure conductive model greatly underestimates the effective 
"conductivity". 

The early theoretical model of Smoot, Horton, aiid Williams^ should 
predict a burning velocity which is much lower than the observed values due 
to its neglect of radiative processes. Yet the model predicts a burning 
velocity which is actually higher than the observed value (see Fig. 4). It 
also predicts a maximum temperature which is higher than the observed value. 

For 10 y particles at concentrations less that 500 mg/& and the lar- 
ger- 33 y particles at all concentrations studied in the experiments (Fig. 4), 


the optical thickness of the flame front is moderate or small. To the 
extent that F is related to the gradient of the mean intensity, 

AJ 


F 'v 


Ax ’ 


the maximum F will remain of the same order of magnitude regardless of the 
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opacity of the flame front. However, their distributions will be different, 

as shown schematically in Fig. 5. For smaller flame front opacity, the 

radiation flux remains important while dF/dx becomes smaller within the 

2 2 

reaction zone in comparison with (k)d T/dx . One expects, therefore, the 
still considerable F in the preheat zone to raise the temperature of the 
dusty medium. Within the reaction zone, the overall exothermic reaction 
continues to increase the temperature through conduction processes. This 
leads to high maximum temperature in the reaction zone. One may confirm 
this conclusion from another line of reasoning. In the preheat zone and the 
post-reaction zone where conduction and heat release are negligible, equa- 
tion 1 reduces to, 


w 


dT 

dx 


iTc S ^ 


or 


T -F . 


With the opacity of the reaction zone being moderate or small, F in the 
preheat zone and the post-reaction zone remains significant (.see Fig. 5). 
Temperatures immediately in front of the reaction zone and behind the reac- 
tion zone will be raised by the combined conduction-radiation process. 

It is therefore expected, that a mere pronounced peak of the maximum 
temperature develops in dust flame when tire opacity of the flame front 
decreases from large to moderate value. The burning velocity should increase 

markedly with decreasing concentration of 10 y particles. This is also in- 
dicated by the observed burning velocity showing a marked increase and be- 
coming greater than the theoretical prediction of a pure conductive model. 

For 33 y particles, the observed burning velocity remains high throughout 
and at 1200 mg/i, it is actually higher than that of 10 y dust. According 



25 


to. the present view, this is due to the pronounced temperature peak which 
remains in the case of 33 y at 1200 mg/Ji. The temperature wanes in the 
case of 10 y particles at 1200 mg/&. 

Finally, we note that the cited experimental data for coal dust flames ^ 
were obtained at normal gravitational conditions. Our experience with (pre- 
liminary) studies of lycopodium dust flames suggests that apparent lean 

limit observations Cat g = 1) do not correspond to those to be obtained under 

24 

gravity-free conditions Also, settling velocities vary with particle size, 
suggesting that (g = 0) vs. (g = 1) observations differ more substantially 
as settling velocities (particle size) increase. It is also to be noted 
that the added sophistication of this conductive-radiative dust flame model 
is particularly appropriate for flames subject to no gravitationally induced 
body forces. This is not the case for the data of reference 6. The theore- 
tical approach discussed here is expected to be directly applicable to data 
derived during Space Shuttle combustion studies. 
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VI . OTHER COMBUSTION PHENOMENA AMD GRAVITATIONAL EFFECTS: 

OSCILLATORY COMBUSTION OF CARBON MONOXIDE . 

2 

In a previous REVIEW we have noted the need for the study of a range 
of other combustion phenomena at reduced gravitational conditions. In the 
following pages, we report on our recent studies in one such area, OSCILLATORY 
COMBUSTION OF CARBON MONOXIDE . In this case we again find that space-based 
combustion experimentation is needed to verify and/or assist in the development 
of existent or needed theory. 

Previous studies of the oscillatory combustion of carbon monoxide have 

2 20 36 

lead to the conclusions that wall effects are of great importance. * 

The physicochemical nature of these wall effects have been the subject of 
substantial discussion, but the essence of these discussions involves the 
transport of heat and mass between the gas phase oxidative processes and the 
reactor walls. In general, it is thought that heterogeneous chemical processes 
at the reactor walls play a key role in defining the oscillatory combustion of 
carbon monoxide. Inasmuch as temperature gradients between the axial regimes 
of such, (long, cylindrical) reactors appear to be small, it is believed that 
heat transfer processes between walls and gas phase are of lesser importance 
than mass diffusion (and heterogeneous kinetic) effects. Accordingly, the 
oscillatory oxidation of carbon monoxide has been thought to be a purely 
kinetic oscillation , as distinct from the well-known thermokinetic oscillations 
found for the oxidation of hydrocarbons. 
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The theory of kinetic oscillations given by Yang and Berlad ’ involves 

( 1 ) 

a 


the following chemical kinetic schematics: 

'^1 

CO + 0, = -> C0 9 + 0 

■ *2 

CO + 0 + M > C0*CC0 9 ) + M 

k ' 

CO* + 0 CO + 0 2 
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CO* + M 
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C0 2 + M 
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The roles of kinetic wall saturation effects and other heterogeneous 
kinetics are implicit in schematics (8), (9), (10), (14), and (15). Unlike 
the (other) rate constants which, describe homogeneous gas phase kinetics, 
kg, kg, h^g, k^, and k^ can be shown, in general, to be functions of tube 
geometry, tube diameter, a diffusion coefficient, and the concentration 
gradients which drives species flux rates between walls and gas phase. Consider, 
then, the diffusive transport of some species, c^, from the gas phase reaction 
volume to an absorbing wall. The diffusion equation, in the absence of natural 
convective effects, can be written 
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( 18 ) 


dc • d 2 C . . r t r 

-t— = D — 7- + c. 
dt .2 1 

: dx 


where, for convenience, the geometric arrangement is taken to be plane parallel 
plates of infinite extents The relation to other geometries will be noted 
later. If the plane parallel plate wall separation is (d) and the reactor 
extends oyer the' range - d/2 < x <^ + d/2 and , for a quasi-steady process, 
equation ( 18 ) can be integrated to yield 


( 19 ) 


. » t 1 

c . 

C i = ~2D~ 


A 2 2 

Gr) - x 


The mean concentration for the reactor is given by 
d/2 / d/2 


( 20 ) 


C i = 


c .dx/ 

1 It 


- d/2 


dx 


-d/2 


.Ml 


Then, assuming that c^ is spatially uniform and that there are no wall re- 
flections, we obtain a relation between the rate of production of the species 
of interest and the other parameters of the problem. 


( 21 ) 


.tit 
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12 Dc. 
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Thus, for example, for the case of schematic number CIO) , 


OH 


v 10 


■p destruction on wall, and 


the corresponding sink term for [OH] implies 

( 22 ) 


V 


10 
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Comparing (21) with (22) shows that 



A more general expression for may be written in the form 



where (Z) is a characteristic geometric parameter having the values 

Z * 12, for plane parallel plates of infinite extent, and 
Z = 32, for cylinders. 

It is significant to note the various assumptions that go into the 
seemingly straightforward schematics which attempt to describe the heterogen- 
eous processes associated with walls. There is the assumption that natural 
convection is absent and that a known quasi-steady concentration (clearly 
nonuniform and unsteady for our case) is operative. Even if one accepts these 

assumptions, the diffusive regime of heterogeneous kinetics depends upon tern- 
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perature, pressure, vessel size and vessel geometry. Gray has shown that 
small nonisothermalities in reactors do lead to natural convective transport 
between walls and reactor volumes. 

Systematic variation of reactor vessel size subject to convection-free 
heat and mass transport Iras not been carried out for the case of oscillatory 
combustion of carbon monoxide. Indeed, if such experiments can be carried 
out under gravity-free conditions, the time and space variations of the various 
species important to carbon monoxide oxidation could be easily followed (e.g., 
in a cylindrical reactor) spectroscopically. The gravity-free case permits 
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true cylindrical symmetry both in the cylindrically uniform heating of the 
reactor vessel as well as the radially symmetric transport of heat and mass. 
Such studies would permit observation of a level of experimental detail 
(Space and time variation of concentrations and fluxes) which have not here- 
tofore been possible. ' 
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VII. LIST OF FIGURES 

Fig. 1. Extinction efficiency, X^ , and ratio of scattering to absorption 
efficiency, X /X , for spherical particles with refractive index 

S cL 

2.00 — 1.20 i (typical of carbon black) as a function of the 
dimensionless particle size (circumf erence/wavelength) . 

Fig. 2. Rosseland mean extinction coefficient for particle size 10 y and 
33 y at 1200°K. 

Fig. 3. (Kp - Kp)/Kp vs. particle size. 

Fig. A. Measured and predicted effects of coal particle size and coal 
■ \ ■ . ■ 

1 ■ concentration on. flame velocity,— — ■ — ; -data correlation, 

theory (monodispersed) . 

Fig. 5. Schematic profiles of the radiation fluxjN , 

op.tically-thick flame; — , optically-thin flame. 
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